Inclusive Ridge Distributions in Heavy-Ion Collisions 
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The formation of ridges induced by semihard scattering in nuclear collisions is included in the 
description of single-particle distributions for both pion and proton at low transverse momenta. The 
ridge component is characterized by an azimuthal dependent factor that is derived in the study of 
the ridge structure in two-particle correlation distributions involving triggers. It is shown that the 
inclusive ridge can reproduce the observed data on V2(pr) if the base component underlying the 
ridge has no azimuthal dependence. A common description of pion and proton spectra is given 
in the recombination model that can smoothly join the low- and intermediate-py regions. All the 
important properties of single-particle distributions in those regions can be satisfactorily described 
in this approach. 

PACS numbers: 25.75.-q, 25.75.Dw 



I. INTRODUCTION 

As the data on single-particle distributions of iden- 
tified hadrons produced in heavy-ion collisions become 
more abundant and precise [lHj| , more demands are put 
on theoretical models to reproduce them. It is gener- 
ally recognized that in Au-Au collisions at -^/saT/v = 200 
GeV at the Relativistic Heavy-Ion Collider (RHIC) the 
low transverse-momentum (pr) region (pr < 2 GeV/c) 
is well described by hydrodynamics |8[ and the high- 
Pt region (pr > 6 GeV/c) by perturbative OCT) [!|, 
both subjects being reviewed recently in Ref. [lfj. In 
the intermediate region (2 < pt < 6 GeV/c) neither 
approaches work very well. What stands out in that re- 
gion are the large baryon/meson ratio and quark- number 
scaling (QNS), which give empirical support to the re- 
combination/coalescence models [lll - tlq . The connec- 
tion between the intermediate- and high-p^ regions is 
smooth, since the dominance of shower-shower recombi- 
nation is equivalent to parton fragmentation. The tran- 
sition across the lower pr boundary at pr ~ 2 GeV/c is 
not so smooth because of the difference in the continuum 
description in hydrodynamics and the parton description 
in hadronization. Our aim in this article is to extend 
our previous considerations [13, EH to the lower-p^ re- 
gion and to describe in a self-consistent way both the pt 
and azimuthal <f> behaviors of pions and protons without 
explicit reliance on hydrodynamics. 

One specific point that motivates our study is related 
to the question of what happens to the initial system 
within 1 fm/c after collision. Semihard partons cre- 
ated within 1 fm from the surface will have already left 
the initial overlap region before thermalization is com- 
plete. There are many of them with parton transverse- 
momentum fcx ~2-3 GeV/c even at RHIC, let alone at 
the Large Hadron Collider (LHC). They are minijets that 
can cause azimuthal anisotropy, not accounted for by con- 
ventional hydrodynamics. It is known that in events trig- 



gered by jets there is a ridge phenomenon in the struc- 
ture of associated particles with narrow Acj) (azimuthal 
angle relative to that of the trigger) and extended Aij 
(pseudorapidity relative to the trigger). Such a struc- 
ture should be present in the inclusive distribution even 
if triggers are not used to select the jet events. When hr 
is low enough so that minijets are copiously produced, 
the corresponding effect on the <j> anisotropy can become 
dominant, rendering the consideration of pressure gra- 
dients along different <f> direction unreliable if semihard 
scatterings are ignored. In this paper we give specific 
attention to the ridge contribution to the single-particle 
distributions in the \ow-pr region. It is in this sense that 
we use the terminology: inclusive ridge distribution. 



Another area of concern is the variation of the pr de- 
pendence as the focus is moved to the low-pr region, 
where pion and proton appear empirically to have differ- 
ent behaviors. In the parton recombination model the 
hadrons should have the same inverse slope as that of 
the coalescing quarks if the hadrons are formed by re- 
combination of the thermal partons, but because of the 
difference in the meson and baryon wave functions, the 
net pt distributions turn out to be different. This line 
of analysis takes into account the quark degree of free- 
dom just before hadronization, which is overlooked by 
the fluid description of the flow effect. The burden is to 
show that the data on V2(pr) can be reproduced for both 
pion and proton at low pr without the hydro description 
of elliptic flow. That is indeed what we shall show for 
various centralities. 



We confine our consideration in this paper to the 
physics at midrapidity. At larger 77 there are other issues, 
such as large p/w ratio [f| and large A-q distribution of 
triggere d ridge [7|, which have been examined in Refs. 
19, 20], and will not be further considered here. 
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II. SINGLE-PARTICLE DISTRIBUTION WITH 
RIDGE 

We begin with a recapitulation of our description of 
single-particle distribution [ll|, [lj], fl8| . At low pt we 
consider only the recombination of thermal partons, so 
the pion and proton spectra at y = are given by 
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where T{qi) is the thermal distribution of the quark (or 
antiquark) with momentum g.;, and 7Z h is the recombi- 
nation function (RF) for h — ir or p. On the assump- 
tion that collinear quarks make the dominant contribu- 
tion to the coalescence process (so that the integrals are 
one-dimensional for each quark along the direction of the 
hadron), the RFs are 
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where the details of /(gj/p) that depends on the proton 
wave function are given in [Tl| , and need not be repeated 
here. The main point to be made here is that if the 
thermal distribution T(<?i) has the canonical invariant 
form 



T (q) =q^= Cqe^ T , 
dq 



(5) 



then the ^-functions in the RFs require that dN h / 'pxdpr 
has the common exponential factor, exp(— pt/T), for 
both h — 7r and p. The prefactors are different; we simply 
write down the results obtained previously 
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where Af v oc C 2 and Af p oc C 3 , and C has the dimen- 
sion (GeV)^ 1 . Note that the factor p\jmT in the proton 
spectrum (that must be present for dimensional reason) 
causes the pjix ratio to vanish as pt — > on the one 
hand, but to become large, as pr increases, on the other. 
When pt exceeds 2 GeV/c, shower partons become im- 
portant and the above description must be supplemented 
by thermal-shower (TS) recombination that limits the 
increase of the p/n ratio to a maximum of about 1 (llj . 

We restrict our consideration to pt < 2 GeV/c, but 
now broaden it to include (f> dependence. For non-central 



collisions the almond-shaped initial configuration leads 
to <fi anisotropy. The conventional description in terms 
of hydrodynamics relates the momentum anisotropy to 
the variation of pressure gradient at early times upon 
equilibration [21(. The success in obtaining the large 
V2 as observed gives credibility to the approach. We 
adopt an alternative approach and justify our point of 
view on the basis that we can also reproduce the em- 
pirical «2, as we shall show. Furthermore, aside from 
offering a smooth connection with the intermediate pt 
region by the inclusion of TS recombination, our ap- 
proach describes also the effect of semihard scattering 
on the soft sector. The ridge phenomenon that we at- 
tribute to that effect can be with trigger [2(1 H^, [23| or 
without trigger 0, E3, El, H3- Although data on the 
ridge structure must necessarily make use of triggers in 
order to distinguish it from background [25-28], inclusive 
distribution must include ridges along with background. 
Thus theoretically a single-particle distribution should 
have a ridge component in the soft sector due to unde- 
tected semihard or hard partons. That component has (j) 
dependence that can be calculated from geometrical con- 
sideration [l8| , and has been shown to be consistent with 
the dependence of the ridge yield in two-particle correla- 
tion on the trigger angle <p s relative to the reaction plane 
27]. 

Let us use Pi(pt, 4>, b) to denote the single-particle dis- 
tribution of hadron h produced at mid-rapidity in heavy- 
ion collision at impact parameter b, i.e. 
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where -/Vp art is the number of participants related to b 
in a known way through Glauber description of nuclear 
collision (29[. At low let p\ be separated into two 
components 



p\{pT, ^ b) = B h ( PT , b) + R h (pt , <t>, b), 



(9) 



where B h (pT, b) is referred to as Base, not to be confused 
with the bulk that is usually determined in hydrodynam- 
ics; this is a change from earlier nomenclature [l8j |. where 
the use of "bulk" did lead to some misunderstanding. 
Our emphasis here is that B(px,b) is independent of <p. 
In our approach we regard the semihard partons created 
near the surface, and directed outward, give rise to all 
the (f> dependence of the medium before equilibrium is 
established; the recoil partons being directed inward are 
absorbed and randomized. The component expressed by 
R h (pT, (f>,b) is referred to as ridge on the basis of its cj> 
dependence discussed below. The B h (pT,b) component 
consists of all the soft and semihard partons that are 
farther away from the surface and are unable to lead to 
hadrons with distinctive dependence. Thus the separa- 
tion between B h (pT,b) and R h (pT,<f>,b) relies primarily 
on the 4> dependence that the ridge component possesses. 

In Ref. [18] we have given an extended derivation of 
what that <j) dependence is. It is embodied in S{4>, b) that 



3 



is the segment of the surface through which a semihard 
parton can be emitted to contribute to a ridge particle at 
0. From the geometry of the initial ellipse (with width 
w and height h that depend on 6) and from the angular 
constraint between the semihard parton and ridge parti- 
cle prescribed by a Gaussian width a determined earlier 
in treating the ridge formation for nuclear density not 
too low I23|, it is found that 



where 



S(4>,b) = h[E(0 2 ,a)-E(O u a)], 



(10) 



where E(8i,a) is the elliptic integral of the second kind 
with a = 1 — w 2 /h 2 and 



tan 1 | — tan < 

w 



a, 2 = + o-(ll) 



for < 7r/2, and an analytic continuation of it for 02 > 
it/2. Thus S((f>,b) is completely calculable for any given 
b, and R h (pT, 0, b) is proportional to it. 

We can now rewrite Eq. Q unambiguously as 
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and R h (p T ,b) is a similar average of R h (pt, 0, b). Ac- 
cording to Eqs. © and ([7]) the inclusive distributions 
Pi(PT,b) should share the common exponential factor 
exp(— pt/T), for h = n and p, as for quarks. That does 
not take into consideration the enhancement of pions at 
very small px due to resonance decay. We account for it 
by a phenomenological term u(pr, b), and write 



pU PT ,b) = NAb)[l + u{p Tl b)]e~ PT/T , (14) 
p\{p T ,b) = M P {b)—z- VTlT , (15) 



where the resonance effect on the proton is neglected be- 
cause of baryon-number conservation. These expressions 
are for the left-hand side of Eq. (fl2|) after averaging. 
The base term B h (pT, b) on the right side is the soft com- 
ponent without the contribution from semihard scatter- 
ing near the surface and should have the same common 
structure as in Eqs. © and ([7]) due to thermal parton re- 
combination, except that the inverse slope is lower with- 
out the enhancement by the energy loss from the semi- 
hard partons. We can therefore write 

B"(p T> b) = K(b)[l+u(p T ,b)]e- pT/TB , (16) 

(17) 



BP(p T ,b) = M P {b)^e~^/ T \ 
rriT 

where Tb < T. It then follows that 
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AT = T — Tb- (21) 



There are two undetermined inverse-slopes: Tb and T, 
common for both 7r and p. They are for single-particle 
inclusive distributions, so only T is directly observable. 
We postpone phenomenology to a later section. In ridge 
analysis using triggered events for two-particle correla- 
tion the two corresponding inverse slopes are separately 
measured [26| . Here, however, we are dealing with single- 
particle distributions. The difference between Tb and T 
has to do with ridges and their effect on the (f> distri- 
bution. Thus we expect AT to be related to azimuthal 
asymmetry, a topic we next turn to. 



III. QUADRUPOLE MOMENTS OF <j> 
ASYMMETRY 

This topic is usually referred to as elliptic flow, a ter- 
minology that is rooted in hydrodynamics. Since we have 
not used hydro in the previous section, it is more appro- 
priate to use the unbiased language initiated in Ref. [30j , 
and call it azimuthal quadrupole. It is the familiar V2 
that is defined by 



J 27r d(j)COs2(j)p l l(p T ,(j},b) 



^(pt,6) = <cos2^ 1 = 



/o' r d#i(PT,0, b) 
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Using Eqs. (H]) - fl5T) yields 
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V2[PT '' ' B(pT,b) + R(p T ,b) 

(cos 20) s 



Z- l {p T ) + l 



(23) 



where 
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(cos 20) s = I d<f> cos 205(0, 6), (24) 
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These equations are remarkable in that the 6 dependence 
resides entirely in Eq. ([24]) and the px dependence en- 
tirely in Eq. (|25[) : furthermore, there is no explicit de- 
pendence on the hadron type nor the resonance term 
represented by u(j>T,b). As we have noted at the end 
of the preceding section, T can be determined by the 
Pt spectra, but Tb is not directly observable. However, 
the quadrupole is measurable, so it can constrain T and 
therefore Tb- In short, the two parameters T and Tb can 
be fixed by fitting the data on Pi(pt, b) and v^Ipt, 6). 

Without using a model to describe the evolution of the 
dense medium, it is clear that we cannot predict the val- 
ues of T and Tb- However, our aim is to discover how far 
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one can go without using such a model. Neither T nor Tb 
depend on 0. Yet non-trivial v^ipTib) can be obtained 
because of the presence of the ridge term in Eq. (JT2|) . 
If phenomenology turns out to support this interpreta- 
tion of azimuthal asymmetry, as we shall do in the next 
section, then the ridges induced by undetected semihard 
partons play a more important role in giving rise to the <f> 
dependence in inclusive single-particle distribution than 
hydro expansion that is based on assuming equilibration 
to be completely at a later time without semihard scat- 
tering. 

From Eqs. PH|) and (IM1) we can calculate (cos 20) s 
and obtain its dependence on b. For the initial elliptical 
configuration the width and height are 



w = 1 - 6/2, h = (1 - 6 2 /4) 



1/2 



(26) 



where all lengths are in units of the nuclear radius Ra- 
Setting the Gaussian width a between the azimuthal an- 
gle 4>i of the semihard parton and 02 of the ridge particle 
to be (7 — 0.33 . we determine (cos 20) s as shown in 
Fig. 1(a). 

According to Eq. (|2"3")l (cos 20) s contains all the b de- 
pendence of v^pxjb) for any px in the soft region. To 
check how realistic that is phenomenologically, we show 
first in Fig. 1(b) the data on v\ (pr, N part ) for three pr 
values from Ref. 0], but shifted vertically so that they 
agree with the data for pr — 0.975 GeV/c for most of 
large N part . The diamond and square points are slightly 
shifted horizontally to spread out the overlapping points 
for the sake of visual distinguishability. The fact that 
their dependencies on N part are so nearly identical is re- 
markable in itself. The solid line is a reproduction of the 
curve in Fig. 1(a) but plotted in terms of N part , and re- 
duced in normalization by a factor 0.25 to facilitate the 
comparison with the data points. For N part > 100 the 
line agrees with the data on vi very well, thus proving the 
factorizability of px and b dependencies of Eq. (|2"3"f . For 
N P art < 100, corresponding to b/R.A > 1.3 or centrality 
> 40%, there is disagreement which is expected because 
the density is too low in peripheral collisions to justify 
the simple formula in Eq. (I23 p. A density-dependent cor- 
rection is considered in Ref. [18J , but will not be repeated 
here. Our focus in this paper is on the inclusive ridge, so 
we proceed to phenomenology on the basis that the for- 
malism given above is valid for central and mid-central 
collisions at N part > 100. To have a compact analytic 
expression for 5(0, b) as given in Eqs. (fT0|) and (jTTT) to 
summarize the dependence is not only economical, but 
also provides a succinct feature to distinguish the ridge 
from the base components in Eq. (11/ 



IV. PHENOMENOLOGY 

We now determine the parameters in our model 
through phenomenology. A success in fitting all the rele- 
vant data can give support to our approach that empha- 
sizes issues not considered in the standard model 1311. 
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FIG. 1: (Color online) (a) Average of cos 20 weighted by 
S(4>,b) vs impact parameter in units of Ra- (b) Common 
dependence of i>2(Pt,&) on N par t for various pr, shifted ver- 
tically for comparison. The diamond and square points are 
horizontally shifted slightly from the points in circles to aid 
visualization. The solid line is from (cos 20) s shown in (a), 
but rescaled and plotted in terms of N part . The data are from 
Ref. [|. 



Our first task is to determine the inverse slope T that 
is shared by T{q),Pi{pT, b) and p\{pt, b). Since the nor- 
malization factors in Eqs. (JSJ), (|14p and (|15l) have not yet 
been specified, we consider first a particular centrality, 
20-30%, and fit the px dependence of the proton spec- 
trum for px < 2 GeV/c, as shown in Fig. 2, and obtain 



T = 0.283 GeV. 



(27) 



Note that the one-parameter fit (apart from normaliza- 
tion) is very good compared to the data from Ref. [l|. 
It demonstrates that the proton is produced in that px 
range by thermal partons and that the flattening of the 
spectrum at low px is due to the prefactor p T /% arising 
from the proton wave function. 

Having determined T, we next consider the pion spec- 
trum p~i(pT,b). According to Eq. fll"4l it has the same 
exponential factor as does Pi(pT,b), but has also an ad- 
ditional factor [1 + u(pT,b)] due to resonance decay. We 
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FIG. 2: Proton spectrum at y ~ averaged over cA (hence, 
no 1/2-7T factor) at 20-30% centrality. The solid line is a fit 
of the data by Eqs. (|15|l and (|27[) with free adjustment of 
normalization. The data are from Ref. 



show in Fig. 3 the data from PHENIX |l| on the pion 
distribution for 20-30% centrality; the exp(— pt/T) fac- 
tor is shown by the dashed line, the normalization being 
adjusted to fit (and to be discussed later). For pt > 1 
GeV/c they agree very well, demonstrating the validity 
of the common T. For pr < 1 GeV/c there is resonance 
contribution to the pion spectrum which we cannot pre- 
dict. Thus we fit the \ow~pt region by the addition of a 
term exp(— px/T r ), shown by the dash-dotted line, cor- 
responding to T r = 0.174 GeV. The sum depicted by the 
solid line agrees with the data perfectly. The point of this 
exercise is mainly to show that the common exp(— px/T) 
behavior is valid for pion as for proton, but the reality of 
resonance contribution for pr < 1 GeV/c obscures that 
commonality. Converting the resonance term to the form 
given in Eq. (fl4|) we write 



u(p T ,b) =u Q (b)e- pT/T °, 



(28) 



where T = 0.45 GeV and u = 3.416 for 20-30% central- 
ity. We do not regard this u term as a fundamental part 
of our model; we attach the factor [1 + u(pT,b)] to all 
expressions of the pion distributions, as in Eqs. (|16[) and 
([!%)) . Of more significance is the role that T has played 
in the phenomenology, and so far Tb has played no role. 

Tb is not directly related to any observable spec- 
trum, since it describes the pr dependence of the base 
B h (pT,b) that lies under the ridge. The important con- 
cept we advance here is that it is <fi independent, and 
that R (pr,<p, b) carries all the c6 dependence. Thus we 
turn to v% (pt , b) in Eq. ([23]) and examine its pt depen- 
dence for both h = 7r and p. In order to emphasize the 
universality between ir and p, we consider v\ versus the 
transverse kinetic energy Et, for Et < 0.8 GeV, where 

E T (p T ) = m T (p T ) - m h . (29) 

We adopt the ansatz that pr is to be replaced by Et in 
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FIG. 3: Pion spectrum showing e _PT//T by the dashed line, 
and the resonance contribution by the dash-dotted line. The 
sum is in solid line. The data are from Ref. [J. 



Eq. (|25|) so as to account for the mass effect, i.e., 

Z{ PT ) = e^ter)/? _ 1; ( 30 ) 

where T is as given in Eq. (|2ip . In Fig. 4 is shown the 
data from Ref. Q when v% is plotted against Et for 20- 
30% centrality. We fit the data points for both h = ir 
and p by Eqs. (l23l) and (|30|) with the choice 

(31) 



T B = 0.253 ±0.003 GeV, 
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FIG. 4: (Color online) v\ for h = n and p. The shaded region 
corresponds to Tb = 0.253 ± 0.003 GeV. The data are from 
Ref. @|. 

which is represented by the shaded region in Fig. 4. The 
upper boundary of that region is for Tb — 0.25 GeV that 
fits the pion V2 almost perfectly, and the lower boundary 
is for Tb — 0.256 GeV that fits well the proton i>2- It is 
evident that vi is very sensitive to Tb due to the expo- 
nential factor in Eq. (I30[) . yet the data support a common 
value for Tb to within 1-2% deviation for pion and proton 
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production. One cannot expect an accuracy better than 
that in the universality of v\ for h = ir and p. We regard 
this result to be remarkable, since the normalization of 



,h : 



is fixed by Eq. 



without freedom of adjustment. 



Note that we have not used any more parameters besides 
T and Tb to accomplish this, which is a fitting procedure 
not more elaborate than the hydro approach where the 
initial condition and viscosity are adjusted. 

So far we have concentrated on 20-30% centrality 
partly because we want to separate the pt and cj) de- 
pendencies from the issue of centrality dependence, and 
partly because v%(pr, b) is large at 20-30% centrality for 
low pt- To extend our consideration to other centrali- 
ties, we fix T and Tb at the values obtained in Eqs. (|27|) 
and (|3"T|) so that Z(pr) is no longer adjustable. The cen- 
trality dependence of " is then examined using Eq. (1231) . 
Figure 5 shows the results for different centrality bins 
for both h = ir and p. The shaded regions due to the 
uncertainty in Eq. pip become narrower in more cen- 
tral collisions. The agreement with data from STAR Q 
is evidently very good. Since there has been no more 
adjustment of free parameters to achieve that, we find 
substantial support from Fig. 5 for our view that the 
<j) dependence arises entirely from the ridge component 
in the inclusive distribution. This raises serious question 
on whether viscous hydrodynamics is the only acceptable 
description of heavy-ion collisions, if the reproduction of 

is the primary criterion for the success of a model. 
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FIG. 5: (Color online) Same as in Fig. 4 for four centrality 
bins. 

It is possible to further improve the agreement between 
the values of v\ for pion and proton in Fig. 5 if those 
figures are replotted in accordance to the idea of quark 
number scaling (QNS), i.e., v^/nt vs Ex/rifi, where rih is 
the number of constituent quarks in hadron h [ljl HH . As 
we have considered QNS and its breaking in the recombi- 
nation model before 17], we do not revisit that problem 
here, especially since our main goal to use vi to constrain 
Tb has already been accomplished. 



Having obtained the correct centrality dependence of 
that is calculable, we now consider the centrality de- 
pendence of the inclusive spectra Pi{pt, b). We note that 
the unknown normalization factors Afir(b) and Af p (b) in 
Eqs. (fT4|) and (fT5|) never enter into the calculation of 
because of cancellation, but for p\{pT,b) they must be 
reckoned with. As remarked after Eq. ([7]), A/^-(6) and 
N p {b) are proportional to C 2 and C 3 , respectively, due 
to qq and qqq recombination. The magnitude C of the 
thermal partons depends on b in a way that cannot be 
reliably calculated. By phenomenology on the pion spec- 
trum it was previously estimated for pt > 1.2 GeV/c 
18], but that is inadequate for our purpose here; more- 
over, Af v (b) and Af p (b) have different statistical factors 
that can depend on b because of resonances. We give 
here direct parametrizations of the normalization factors 



in terms of N, 



part 



Np{N p art) 

uo(N part ) 



= 0.516^, 
= 0.149^1 
= 2.8 + 0.003AT part . 



(32) 
(33) 
(34) 



The parameters are determined by fitting the centrality 
dependence to be shown, but the essence of our predic- 
tion is in px and <f> dependencies that are not adjustable. 
Using the above in Eqs. (fT4")l and (TT5"j) we obtain the 
curves in Fig. 6 (a) pion and (b) proton for three cen- 
trality bins. They agree with the data from PHENIX [l[ 
very well over a wide range of low px- In all those curves 
T is kept fixed at 0.283 GeV, thus reaffirming our point 
that both pions and protons are produced by the same 
set of thermal partons despite the apparent differences in 
the shapes of their pt dependencies. 



V. INCLUSIVE RIDGE DISTRIBUTION 

It is now opportune for us to revisit the two-component 
description of the single-particle distribution and focus 
on the ridge component, in particular. As stated explic- 
itly in Eq. (IT2l . the <f> dependence separates the B h (px, b) 
and R h (pT,b) components, the former being described 
by Eqs. CGI and {T7|), the latter by Eqs. (JI5J) and (E 
Upon averaging over 0, we have 



p1(p T ,b)=B h ( PT ,b) + R h (p T ,b). 



(35) 



Since the exponential factors are the same for h = n and 
p, let us consider only the pion distribution specifically. 
In Fig. 7 we show B and R components by dashed and 
dash-dotted lines, respectively, for (a) 0-5% and (b) 20- 
30%. It is in those figures that we exhibit the basic dif- 
ference between our description of inclusive spectra and 
those of others. Inclusive ridge represented by R is al- 
ways present in the single-particle distribution whether 
or not an experiment chooses to do correlation measure- 
ment to examine the ridge. Semihard scattering is un- 
avoidable in any nuclear collisions at high energy. Its 
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FIG. 6: Inclusive spectra at three centralities for (a) pion and 
(b) proton. The data are from Ref. 



effect on soft partons is therefore also unavoidable. We 
quantify the effect by the R component which is deter- 
mined by the azimuthal anisotropy that is well repro- 
duced in Fig. 5. Here in Fig. 7 we see it rising above the 
^-independent base B component when pt is higher than 
1.4 GeV/c. It is a consequence of the recombination of 
enhanced thermal partons. For px > 3 GeV/c in addition 
to the inclusive ridge the jet component of the semihard 
partons themselves manifests in the spectra in the form 
of thermal-shower recombination that characterizes the 
intermediate-pT region. Thus we have a smooth transi- 
tion from low- to intermediate-pr regions by recognizing 
the importance of the inclusive ridge component. 

It is observed that the dash-dotted lines in Fig. 7 are 
not exactly straight because the ridge component is not 
exponential in pr . However, for pt > 1 GeV/c, R n (pT,b) 
can be well approximated by pure exponential. In Fig. 8 
we show by the solid line the pt dependence of Rq(pt), 
defined in Eq. ([20)1 ; it is the part of the ridge distributions 
R h (pT, b) in Eqs. (fl"8f and (fT9j) that is common for h = it 
and p and is independent of b. From the values of T and 
Tb that we now know, we have T = 2.39 GeV. The (red) 
dashed line is a straight-line approximation of the solid 



FIG. 7: Inclusive distributions for pion showing the base (B) 
component by dashed line and ridge (R) component by dash- 
dotted line for (a) 0-5% and (b) 20-30%. The solid line is 
their sum. The data are from Ref. p[. 



curve for pt > 1 GeV/c by 



Ro(pt) ~ Roe 



-pt/t' 



T' = 0.326 GeV. 



(36) 



Thus the ridge distribution is harder than the inclusive 
distribution characterized by T = 0.283 GeV. This is a 
property that is known from triggered ridges (26|, but 
now it is for untriggered inclusive ridge. 

The enhancement of T' over T is an important point 
to note. Physically, it means that the ridge is a conse- 
quence of the passage of semihard partons through the 
medium, whose energy losses enhance the thermal par- 
tons in the vicinities of the trajectories. We know that 
the enhancement factor is Z(j>t), which has the neces- 
sary pt dependence to render ufCPTjfr) to be in good 
agreement with the quadrupole data. In particular, the 
property that Z{pr) — > as pr — > is essential to guar- 
antee that v^ipTjb) — > in the same limit. The effect 
of Z(px) at larger pt is to increase Tb to T". Although 
Z(pt) increases exponentially, its net effect on R (pt) 
is suppressed by e~ PT ^ T]B . The effective inverse slope T" 
for pr > 1 GeV/c is larger than T of the inclusive by 
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0.5 1 1.5 2 



p T (GeV/c) 

FIG. 8: (Color online) The pr dependence of Ro(pt) defined 
in Eq. (|20f) . represented by the solid line. The (red) dashed 
line is a straight-line approximation for pr > 1 GeV/c, ex- 
pressed by Eq. (|36}. 



43 McV, roughly the same as what Putschke reported on 
the first discovery of ridge, where the triggered ridg e has 
T" larger than that of the inclusive by 45 MeV [331 ] . 

There are, however, subtle differences between trig- 
gered and untriggered ridges. Experimentally, it is nec- 
essary to do correlation measurements to learn about the 
properties of the ridge, which is extracted by a subtrac- 
tion scheme. The inverse slope T" can be compared to 
the inclusive T of the background. For single-particle dis- 
tribution the only measurable inverse slope is T for the 
inclusive spectrum. Theoretically, we assert that ridges 
do not disappear just because triggers are not used. The 
inverse slope T' for Rq(pt) cannot be measured directly. 
It is larger than both T and Tb because Ro(pt) is the 
difference between the two exponentials for pi and B, 
represented by the middle term in Eq. (|20[) . which van- 
ishes as pr — > 0. A physically more sensible way to com- 
pare the various inverse slopes is to recognize that T' is 
significantly larger than Tg because of the enhancement 
effect due to semihard scattering, and that T is the ef- 
fective slope of the inclusive distribution, B + R, that is 
measurable and is between Tb and T' . 

Although our concern in this paper has been restricted 
to the midrapidity region, the physics of inclusive ridge 
can be extended to non- vanishing pseudo-rapidity rj. In 
Ref. [13] a phenomenological relationship is found be- 
tween the triggered ridge distribution in A77 and the in- 
clusive distribution in r\ with the implication that there 
is no long-range longitudinal correlation. However, there 
can be transverse correlation due to transverse broaden- 
ing of forward (or backward) soft partons as they move 
through the conical vicinity of the semihard partons. The 
enhancement of the thermal partons due to energy loss 
is just as we have described in this paper. Indeed, the 
term representing the enhanced px distribution in Ref. 
[20} is essentially identical to that expressed in Eq. (I2TH) . 



Similar consideration has also been used in the explana- 
tion of the ridge structure found at LHC [H], HH] . In this 
paper we have presented the most detailed quantitative 
analysis of the RHIC data in the formalism of inclusive 
ridge that sets the foundation for the ridges at \Ar)\ > 0. 



VI. CONCLUSION 

Our study of inclusive ridge distributions has consoli- 
dated earlier exploratory work with firm phenomenolog- 
ical support, and therefore succeeded in extending the 
hadronization formalism from intcrmediate-pr region to 
below 2 GeV/c, exposing thereby an aspect of physics 
that has not been included in other approaches. The ef- 
fect of semihard scattering on soft partons is accounted 
for by the ridge component whose azimuthal behavior is 
totally characterized by S(<fi, b); it is a calculable quantity 
bounding the surface segment through which semihard 
partons can contribute to the formation of a ridge particle 
at (f>. In an earlier paper [l8| we showed the connection 
between S(<fr, b) and the dependence of the triggered ridge 
yield on the <j> s of the trigger angle relative to the reac- 
tion plane. Now, we have exhibited the central role that 
S(4>, b) plays in determining the azimuthal quadrupole 
vtj:(pT,b) of inclusive distributions. Thus the inclusive 
ridge distribution that we have advanced in this series of 
work serves as a bridge between the single-particle distri- 
bution and the two-particle correlation. Since semihard 
partons are copiously produced before thermalization is 
complete, it is an aspect of physics that should not be ig- 
nored. The success in fitting v^ipr, b) for all central and 
mid-central collisions and for both tt and p by one pa- 
rameter Tb therefore leads to claim of relevance as much 
as viscous hydro does. 

Another attribute of our approach is to unify the pro- 
duction of pions and protons in one hadronization scheme 
based on the recombination of enhanced thermal partons 
so that their spectra have the same inverse slope T de- 
spite apparent differences in the low-py data. That same 
scheme when extended to pr ~ 3 GeV/c explains readily 
the observed large p/ir ratio. Since the ridge components 
in our formalism are the same for it and p, we can predict 
that the p/ir ratio is also large in the triggered ridge as 
in the inclusive. 

The property about ridge that most investigators are 
concerned about is the large Arj range found in correla- 
tion experiments. That is an aspect of the problem that 
has been addressed in Ref. [2(|. Our focus in this paper 
is on the hidden aspect of the ridge that is not easily de- 
tected, but is pervasive because it is in the inclusive dis- 
tribution. Phenomenological success found in this paper 
puts the idea on solid footing. If the concept of inclusive 
ridge is important at RHIC, then its relevance at LHC 
will be predominant. Since the structure of v^ipr, b) ex- 
pressed in terms of S((j>, b) in Eq. (|23p is independent of 
initial density, viscosity, or even collision energy, except 
T, we would expect v-i measured at LHC to be essentially 



9 



similar to what is shown in Fig. 5 for both tt and p. A 
preliminary look at the data from ALICE [36] leads us to 
believe that such an expectation may not be unrealistic. 
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